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We present a new Liouville-integrable natural Hamiltonian system 
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1 Introduction 

A Hamiltonian system is called natural if its Hamiltonian is a sum of a 

Cs| ' positive-definite kinetic energy and a potential. Natural Hamiltonian systems 

^ ■ on (cotangent bundles of) closed surfaces admitting integrals polynomial in 

Tjj- ! momenta are interesting for a number of reasons: 
O 

^ ■ 1. They are classical: The formulation of the problem is due at least to 

c3 ■ Darboux jTj. 



2. If a natural Hamiltonian system admits a real- analytic integral, then 
it admits an integral polynomial in momenta 



3. In particular, all known natural integrable Hamiltonian systems on 
surfaces have integrals polynomial in momenta. 
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4. The existence of the integral polynomial in momenta of degree one 
or two has a very clear geometric background: the existence of an 
integral of degree one implies the existence of a one-parametric family 
of symmetries of the system. The existence of an integral of degree two 
implies the existence of so-called separating variables. 

Natural Hamiltonian systems on closed surfaces that admit a nontrivial 
integral which is polynomial in momenta of degree one or two are completely 
understood: there exists a complete description and classification, see, for 
example, j3]. Up to now, only one family of natural Hamiltonian systems on 
a closed surface admitting an integral of degree three and admitting no (non- 
trivial) integral of degree one or two is explicitly known, namely the family 
of system obtained from the Goryachev-Chaplygin top. We now present a 
second, distinct, family of natural Hamiltonian systems with a cubic integral. 

1.1 The result 

Consider the sphere S"^ C M.^ of radius 1 with the spherical coordinates 
{x,y,z) = {— sin 6 cos (f), — sin 6 sin (f), cos 6) and the following two functions 
H (the Hamiltonian) and F (the second integral) on the cotangent bundle of 
the sphere without poles z = ±1. Let A, c, s G M be parameters with s > 1 
and define 

W{z) = z + s, P{z) = 3z'^ + 4:sz + l, Q{z) =3z'^ + 2SZ-1 

and 

G(.) «^' 



{2W{z)y' 
Then, the function H is given hj H := K + V , where 

TA ^ sin6' , c 

V := A , = cos + 



Vw^(cos^) w{cosey 

and F is defined as 

F := 2Hp^-pl + Acos(0)-= |^\ p^ + 2Asinct>^W{cose)pe . 

v/iy(cos6')sm6' 

Proposition: The functions H and F can he analytically continued to the 
cotangent bundle of the whole sphere. The continuation is also polynomial in 
momenta (of degree 2 for H and of degree 3 for F). 



We will denote the continuations of if, K, F by the same letters if, K, F; 
in particular in the theorem below we mean the continued functions defined 
on the cotangent bundle of the whole sphere. 

Theorem: The following statements hold: 

1. The functions H and F commute with respect to the standard Poisson 
bracket on T*S'^ and are functionally independent 

2. The kinetic energy K is positive definite 

3. If A 7^ 0, the Hamiltonian H does not admit a (smooth) nontrivial 
integral which is polynomial in velocities of degree less than three and 
which is linearly independent of H. 

In other words, we found a new 1) Liouville integrable, 2) natural Hamil- 
tonian system on the sphere with 3) an integral cubic in momenta. 

1.2 Why the system is new 

The previously known family comes from the Goryachev-Chaplygin top using 
symplectic reduction, see P] for the details. It is a natural Hamiltonian 
system on the sphere. In the standard spherical coordinates its Hamiltonian 
iii is Ki + Vi, where the kinetic energy Ki and the potential Vi are 

1 /C0S2^ \ 2 1 2 

Vi := Ai sin 6* sin , 

see, for example, [0]. The integral Fi for the Hamiltonian iii is given by 

. Ai , , / X ^1 sin0(3 (cos^)^-2) 

Fi := iiip,^ - 2p^ - -—pe cos{(j)) cos(e) + — ^— p^. 

z z sm (7 

It is easy to see that the system we found is essentially different from 
the Goryachev-Chaplygin system, i.e. that there exists no diffemorphism D : 
S^ — ^ S^ such that D^K = aKi and D^^V = [3Vi for certain constants a, /3. 

Consider the metrics g, gi corresponding to the kinetic energies K, Ki, 
respectively. They are given by the formulae: 
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By direct calculation, it is easy to see that level lines of the curvature of 
both metric are the lines {6 = const}. More precisely, since the formulae 
above do not depend on 0, the curvature is constant along every such line. 
It is easy to see that the curvature is not constant. 

Then, the diffeomorphism D, if it exists, must take the lines {6 = const} 
of the first system to the lines {9 = const i} of the second system. Since 
the lines {0 = const} are orthogonal to the lines {(f) = const} in the first 
and in the second metric, the diffeomorphism D must also take the lines 
{(j) = const} of the first system to the lines {0 = const i} of the second 
system. Using that the coordinate is periodic with the period 27r, we 
obtain that the diffeomorphism D must be given by formulae 

9new = 9new{9), (pnew = ±0 + 4>0 ■ 

Then, the pull-back of the metric gi has the form 

dO + 



09 J 4 + cot2^„ 

Then, ^^^^ = ±1, and the diffeomorphism D must be given by the formula 

Onew = 6 or 6new = 71-9, (j)new = ±0 + 00 • 

Clearly, the pull-back of the metric gi with the help of such diffeomorphism 
is not proportional to g. Otherwise G would need to become a constant. But 
then P would have a double root at z = —s, which is impossible. Thus, our 
system is essentially different from the reduced Goryachev-Chaplygin top. 

Selivanova ^Tj proved the existence of an additional family of natural 
Hamiltonian systems admitting integrals of degree 3 in momenta. This family 
is not explicit. Instead, a nonlinear differential equation is derived and it is 
proved that certain solutions of this equation allow to construct a geodesic 
flow admitting an integral of degree three and admitting no integral of smaller 
degree. 

At present we do not know whether our system overlaps with the systems 
considered by Selivanova. Our guess is that this is not the case. Here are 
some reasons supporting this claim: 

1. It is not clear whether the system of Selivanova is real-analytic; she 
claims the smoothness only. 

2. Selivanova's family contains (the reduction of) Goryachev-Chaplygin 
system; our family does not. 

3. The ansatz of Selivanova is more restrictive than the form of our system. 



2 Proof of the proposition 

2.1 Analytic proof 

Clearly, the functions V, K, F depend real-analytically on (0, 6) and on the 
corresponding momenta p^,Pe, so the only possible problem could appear 
near the points z = ±1. We consider (x, y) as a local coordinate system near 
these points of the sphere. The coordinates (x, y) are related to coordinates 
{(f), 6) by the following real-analytic change of variables: 



sin 9 = a/x^ + y"^, pg = {xp^ + ypy) ^J [x^ + y"^) ^ - 1 

tan = y/x, p^ = xpy - yp^ . 

In the new local coordinates the kinetic energy and potential are 
K = -{pI+pI- {xpa, + ypyf + {xpy - yp^fG{z)) , 



X 



V{x,y) = —A — -^=^ + , where z = ±\/ 1 



x"^ — y"^ . 



vw{^ ^i^y 

The sign of z depends on whether a neighbourhood of the north or south pole 
is considered. Clearly, K and V are real-analytic. Thus, the Hamiltonian H 
is a real- analytic function on the whole T*S'^. Since we will need it for the 
proof of the second statement of the Theorem, we remark that the kinetic 
energy is positive definite near the poles, where x and y are small. 

Now let us show that the second integral also is a real-analytic function 
on the whole T*S'^. The first two terms of F are real- analytic, since the 
momentum p^ is xpy — yp^ in the new coordinate system, and therefore is 
real- analytic, and we already have proven that H is real- analytic. The sum 
of the third and the fourth terms in the new coordinate system is 




-xy 2 - 3a;^ - 2y' + 2sWl - x^ - y^ 

Px H , Py 



'1 — x"^ — y"^ + s y a/1 — x^ — y^ + s 

and this is clearly real-analytic. 

2.2 Geometric proof 

The second integral F is cubic in the momenta, like the Goryachev-Chaplygin 
top. This top does not appear as a limiting case in our family, but is serves as 
a motivation for the following coordinate transformation, which uses global 
coordinates {x,y,z) on a sphere embedded in R^. The angular momentum 



(Lj., Ly, Lz) where L^ = xpy—ypx = p^p etc. with cychc permutations, satisfies 
the usual Poisson structure on the sphere with non-vanishing brackets 

{Lx, Lyj = Lz, {x, Ly\ = z, {z, Ly\ = —X 

and cychc permutations thereof. This bracket has Casimirs x'^ + y"^ + z"^ = 1 
and xLx + yLy + zLz = 0. The global Hamiltonian is H = K + V, where 

K = l{Ll + Ll + {l + G{z))Ll) 

X c 

V{x,y,z) = -A- 



\/W{z) WizY 
and is clearly real-analytic. The integral reads 

F = 2HLz -Ll- A^^ - 2AVw y^^''y " f ^ 



A 
2HLz -Ll + ^={xLz + 2WLx) . 



In the second line the apparent singularity in F is removed by using both 
Casimirs x"^ + y"^ + z"^ = 1 and xL^ + yLy + zL-, = and in addition the 
identity Q = 2zW + z"^ — 1. The integral is cubic in the momenta, and 
clearly is real-analytic. 

Like the Goryachev-Chaplygin top the system is only integrable on the 
level set of the Casimir xL^ + yLy + zL^ = 0. Writing the Hamiltonian 
in this form suggests an interpretation of a spinning top with zero angular 
momentum (the vanishing Casimir), an orientation dependent moment of 
inertia, and a non-standard potential. Note, however, that 1 + G is not 
positive for all z G [—1, 1] for certain values of s, which means that a moment 
of inertia may pass through infinity for such value of s. Nevertheless, the 
kinetic energy in H is positive definite for all s, see Section 13.21 This is 
possible because even when 1 + G < because xL^ + yLy + zL^ = 0. 

Introducing ^ = x + iy and r] = L^ + iLy and writing V{x, y, z) = xU{z) 
the equations of motion are 

-t^ = -Lz{l + G)^ + zrj (3) 

-tV = -\G'Lli - LzGv + zU- hj' {e + il) ■ (4) 

Here z^ = 1 — ^^ and —2zLz = ^f] + rj^, or the equations are completed with 
z = xLy — yLx and L^ = —yU. 



3 Proof of the Theorem 

3.1 Liouville integrable 

It is tedious but easy to check that the functions H and F commute, since 
they are given by exphcit formulae. The calculation goes as follows. The 
canonical bracket is a quadratic polynomial in the momenta. The term inde- 
pendent of the momenta vanishes because Q = 2zW — 1 + z^. The coefficient 
of pep<j) vanishes because P = 4zW + 1 — z"^. The nontrivial fact is that the 
coefficients of p| which contains P, Q and derivatives thereof vanishes. 

H and F are independent because H has terms quadratic in momenta 
(and not a square of a linear function of momenta) and F has terms that are 
cubic in momenta (and not a square of a linear function of momenta). So 
H and F have vanishing Poisson bracket and are functionally independent, 
hence the system is Liouville integrable. 

3.2 Natural 

As we already have remarked in Section 12.11 the kinetic energy K is positive 
near the poles of the sphere. Thus, it is sufficient to show that the metric 
(P) is positive definite at regular points of the spherical coordinate systems, 
i.e. we have to show that C := 1/ sin^ 9 + G{cos9) is positive for < 9 < n. 
Substituting z instead of cos 9, we obtain 

16z^ + 12sz + 4s2 - 3z^ - Asz^ - 1 

^^^ ~ 4 {i~z^)(z + sy ■ 

Since the denominator {1 — z"^) {z + s) is always positive for —l<z = 
cos (6*) < 1, we need to prove that the numerator 

R(z) := Iz^ + ^sz + s^- -z^ -sz""-- 

is positive. Its value at z = —1 is {s — 1)^ and hence positive; its derivative 
is 

Since the derivative is also positive for —l<z = cos{9) < 1, R{z) is positive 
for —l<z = cos{9) < 1. Hence the metric (PJ is positive definite. 

3.3 Cubic 

Let us prove that, ioi A ^ 0, our system does not admit an integral linear in 
momenta, and that every integral quadratic in momenta is proportional to 
the Hamiltonian. 



Suppose the function P2 + Pi + Pq, where every Pi is a homogeneous 
polynomial in momenta of degree i, is an integral ioT H = K + V. Then, the 
Poisson bracket {K + V, P2 + Pi + Pq} must vanish. We have 

= {K + V,P2 + Pi + Po} 
= {K, P2} + {K, Pi} + {K, Po} + {V, P2} + {V, Pi} + {V, Po} . (5) 

Since the Poisson bracket of two functions -Fi, -F2 '■ T*M'^ — > M is given by 

dx dpx dy dpy dpx dx dpy dy 

where x,y are local coordinates on M^ and Px,Py are the corresponding mo- 
menta, every term on the right-hand side of (jSj) is polynomial and homoge- 
neous in momenta: {K, P2} of degree 3; {K, Pi} of degree 2; {K, Pq}+{V, P2} 
of degree 1 and {V,Pi} of degree 0. The Poisson bracket {y^,Po} evidently 
vanishes. For (jSJ to vanish all the other degree 0, 1, 2, and 3 polynomials 
must also vanish identically. 

Let us prove that the linear term Pi must be zero. Vanishing of {K, Pi} 
implies that Pi is a linear integral of the geodesic flow of the metric given 
by K. In the spherical coordinates 0, 6 on the sphere the metric is given by 
(P) . We have already remarked in Section 11.21 that the level curves of the 
curvature are the curves {6 = const}. 

Suppose the integral Pi has the form Vi{9,(f))p0 + V2{9,(f))p^. Then, by 
Noether's Theorem, the vector field v := (fi,f2) is a Killing vector field, i.e. 
its flow preserves the metric. Then, the flow of v preserves the curvature, 
which implies that the component Vi must be zero. 

Now, using that {V, Pi} vanishes, we have that {V,f2P<^} is zero, which 
in the case f 2 ^ implies ^ = 0, which is definitely not the case. Thus, 
t>2 = 0, and Pi = 0. 

Now let us show that P2 + Pq is proportional to the Hamiltonian H. 
Since {K, P2} = 0, P2 is an integral of the geodesic flow for the metric 
dH). As we already mentioned, for geodesic flows on closed surfaces the 
integrals quadratic in velocities are completely understood. In particular, 
if the curvature of the surface is not constant, then the space of quadratic 
integrals is at most two-dimensional fl31 E] • 

The curvature of the metric (0) is not constant, and the following two 
linearly independent functions are integrals of its geodesic flow: the first is 
K itself, and the second is p'i. Then, for certain constants /?, a, we have 
P2 = ap"^ + j3K. Our goal is to show that P2 + Pq is proportional to H. 
/.From {K, Pq} + {V, P2} = we find {K, Pq - i5V} + a{V,pl} = 0. Hence, 



d{Po-m , 1 djPo - (3V) dV 

oO sm 6 + G{cos6) d<p ocp 



Therefore ^^\f^^ = 0, and then ^^\f^^ = 2a{sm^ 6 + G(cos^))|^. By 



ae 
cross differentiation we find 



d({smH + G{cose))'- 

'09 ""' 



which is possible only if a = or A = 0. If a = then Pq — (3V = const, 
and so after absorbing the constant in the Hamiltonian we have (3 = 1 and 
the integral P2 + Pi + Pq is proportional to H. This proves the Theorem. 

4 Geodesic flows with cubic integrals 

A geodesic flow is a natural Hamiltonian system whose potential energy is 
identically zero. 

According to Maupertuis's principle, an integrable natural Hamiltonian 
system immediately gives a family of integrable geodesic flow, see [3^ for 
details. If the integral of the system is polynomial in momenta, the integrals 
of the geodesic flows are also polynomial of the same degree. 

In our case, the Hamiltonian of the geodesic flow is given by 



^geod 



2{h-V) 



Here /i is a constant such that h is larger than the maximum of the potential 
V{x) on the sphere S"^. The integral of the third degree is 

Fgeod =2Hp^ - 2Vp^ + 4 + Hgeod (2Vp^ 

0(cos9) \ (^) 

+ A cos(0) , J/ J , -p^ + 2 A sin (Py^W {cos 9)pe ^ 

The metric corresponding to Hgeod is 



"•^geod 



{h - V) (d9' + , ,'^'^\ -^ . (7) 

^ ' V l/sin2e + G'(cos^)y ^ ' 



The investigation of geodesic flows that admit an integral polynomial in 
momenta is a very classical subject. If the degree of integrals is one or two, 
they were well-understood already in 19th century, see for example |7, . There 
are a lot of local examples of geodesic flows admitting integral polynomial 
in velocities of degree 3, see the survey PJJ for details, and only very few 
global, i.e. on closed or complete surfaces. 
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Let us recall the known results about geodesic flows on closed surfaces 
admitting an integral polynomial in momenta. For a more detailed review 
see jS]. 

First of all, in view of results of Kolokoltsov jT3], a geodesic flow of a 
surfaces of genus greater than two cannot admit a nontrivial integral poly- 
nomial in momenta. Then, an orientable surfaces admitting such geodesic 
flows must be the sphere or the torus. 

We collect the main results about existence and classification of metrics on 
the torus and on the sphere whose geodesic flow admits a nontrivial integral 
polynomial in momenta in the following table: 





Sphere 5*^ 


Torus T2 


Degree 1 


All is known 


All is known 


Degree 2 


All is known 


All is known 


Degree 3 


Series of examples 


Partial negative results 


Degree 4 


Series of examples 


Partial negative results 


Degree > 5 


Nothing is known 


Nothing is known 



The following notation is used in the table. "Degree" means the smallest 
degree of an integral polynomial in velocities. "All is known" means that 
there exists an effective description and classification (can be found in |2|). 

There exists only one explicit "Series of examples" for degree three. It 
comes from the Goryachev-Chaplygin case of rigid body motion (by applying 
symplectic reduction and Maupertuis's principle), see 0. The most valuable 
"partial negative results" are due to Byalyi JO^ and Denisova and Kozlov [H]. 
They proved that if a natural Hamiltonian system on the torus whose kinetic 
energy is given by a flat metric admits an integral polynomial in momenta 
of degree three, it admits an integral linear in momenta. 

Thus, our system gives one more series of examples of integrable geodesic 
flows on the sphere whose integral is polynomial in momenta of degree 3. 
Below, in Sections 14.1114. 21 we show that these geodesic flows do not admit an 
integral which is polynomial in momenta of degree less than 3, and that these 
examples are different from examples coming from Goryachev-Chaplygin. 
Note, that this does not follow directly from the Theorem. 

Kiyohara fI3] proved the existence of one more series of examples of inte- 
grable geodesic flows on the sphere whose integral is polynomial in momenta 
of degree 3. Although his examples are not explicit (similar to Selivanova 
fn\ he writes a nonlinear differential equation and claims that its solutions 
give integrable geodesic flows), it is easy to see that his examples are different 
from ours, since they have regions of constant curvature and therefore can 
not be real- analytic. 
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4.1 Cubic 

We need to prove that the Hamihonian Hgf,od has no integral which is poly- 
nomial in velocities of degree 1,2. We will prove it assuming that the energy 
level h is sufficiently large. Our proof for arbitrary h is based on the technique 
developed in [HI and is very long. It will be published elsewhere. 

It is easy to see that the integral Fg^od is not the third power of a function 
linear in momenta and is not the product of a function linear in momenta and 
the Hamiltonian. Then, if FaddiUonai ^ is linear or quadratic in momenta 
integral for Hgeod, it must be functionally independent of Fgeod- Then, the 
system is resonant. Then, all geodesies are closed. The theory of metrics all 
whose geodesies are closed is well-developed, see p. It is known [10] that, 
on the two-sphere, geodesies of such metrics have no self-intersections. But 
if h is very big, certain geodesies of our metric must have self-intersections: 
indeed, for huge h, the geodesies of our metric are very close to the geodesic of 
the metric (^, and by direct analysis it is easy to check that certain geodesies 
of (Q) do have self-intersections. 

Thus, for very big h, the Hamiltonian Hg^od admits no integral of degree 
one or two. 

4.2 Why <is^gp^ is new 

The examples coming from Goryachev-Chaplygin are the metrics 



,4^0„^v,)[df-^^^^^]. (8) 



Their Hamiltonians are 



Ho: 



2(/ii - Vi) 
and the integral for the system with the Hamiltonian H2 is 

F2 :=Hip^ - 24 - Vip^ 

M A/ ^1 fA^ rm M ^1 sin0(3(cosgf-2) 

\ 2 2 smB 

Our goal is to show that there is no diffeomorphism D : S"^ -^ S"^ that 
takes the metric (|7j) to the metric (jH)) or a constant multiple thereof. Since 
the potential energy Vi depends linearly on the parameter Ai, it is sufficient 
to show that there exists no diffeomorphism D : S"^ —>■ S"^ that takes the 
metric ((H) to the metric (jH)). 

11 



First let us show that such diffeomorphism D must be given by the formu- 
lae Onew = ^newi^), 4>new = ±0 + 00) where (pQ is a coustaut. It is known, that 
the Goryachev-Chaplygin system is nonresonant on every energy level, see 
for example p!^ where the rotation function is calculated explicitely. Then, 
the geodesic flow of (jH)) is nonresonant as well. Then, the diffeomorphism D 
must take the integral (jH)) to the integral (jHI) or a constant multiple thereof. 

It is known that a Riemannian metric on a (oriented) surface defines a 
complex structure on it. That is, in a neighbourhood of every point we can 
find a coordinate z = x + iy such that the metric looks like X{x, y){dx'^ + dy'^). 
In such complex coordinates, an integral of degree three has the form 

a{z)p^ + b{z)p^p + b{z)pp^ + a{z)p^. (9) 

Here a, b are functions; "bar" means the complex conjugation, p is the com- 
plex momentum corresponding to z. The functions a, b are a-priory not 
assumed to be holomorphic. Note that this form is not restrictive: If © is a 
real- valued function, the coefficients at p^ and p'^p must be complex-conjugate 
to the coefficients at p^ and p'^p, respectively. 

We will use the following result of Kolokoltsov from ^21; in a weaker 
form, this result has been already known to Birkhoff j2]. 

Lemma: Suppose ^ is an integral for the geodesic flow of g and that it is 
not the product of the Hamiltonian and an integral linear in velocities. Then, 



^ dz"^ (10) 



a{z) 
is a meromorphic {3,0) -form. 

Evidently, the form (fTUj) has no zeros. It could have poles. In fact, because 
of Abel's Lemma, it always has poles on the sphere. 

The metrics (I7l8p and the integrals Fgeod, -^2 are given explicitly. Direct 
calculation of the form (fTUj) shows that the form ()10|) for the metric (|7j) and 
for the integral Fgeod is very similar to the form (fTUI) for the metric (jHl) and 
for the integral F2: they have precisely two poles. The poles are located at 
9 = 0,TT. Then, the diffeomorpism D takes the points 6' = 0, vr to the points 
^ = 0,7r. 

It is known that a holomorphic diffeomorphism of the Riemann sphere 
C that preserves the points z = and 2; = 00 is a linear transformation 
z ^— *> az. A holomorphic diffeomorphism that interchange the points is a 
transformation z t^ -. In both cases, the diffeomorphism D has the form 

&new = 0new{d), 4>new = ±0 + 00- 
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Then, the pull-back of the metric (jHI) has the form 
(/ii- Asin^„,^sin(±0 + 0o)) I ( ^1 de' + 



de / 4 + cot2( 



'new I 



If this metric coincides with the metric (|7j), the coefficients of dO^ must co- 
incide. Then, 

de„,,,„V A ■ la \ ■ i^A s A.\\ h ^sin^cos(/) c 

[hi - Ai sva.[enew) sm(±0 + 0o)) = h 



which implies (for simplicity, we assume sign(74) = sign(y4i), the case sign(y4) 
— sign(y4i) can be triated similarly) 

dOnew \^ Asm 6 



^-'■'<*"™) I iT j -^^mTs' <"' 



de„ ^' 



'" I iT j H" - S^WTI 1 ■ (^^' 



Similarly, comparing coefficients of (i0^, we obtain 
Ai sin(6'„e«,) _ Asin(6') 

4 + COt^ Onew ./cOs(^yTs ' 

hi h — c/{cos9 + s) 



4 + cot^ Onew 1/ sin^ e + G{cos 6) ' 
Comparing Equations (J11I13J) and Equations (jl2|14j) . we obtain 

0"new \ . , , 2 n i ' OUnew \ 4 + COt (7^ 

4 + cot Unew and 



(13) 



(14) 



d9 J "^'" \ de J l/sin^^ + G(cos^) 

Thus, we obtain a contradiction, since clearly 1/ sin^ 9 + G{cos 6) is not con- 
stant. Finally, there exists no diffeomorphism D that takes the metric ((Tj) to 
metric (IHI). 



5 What to do next 

We constructed a new natural integrable Hamiltonian system on the sphere. 
It suggests the following ways of further investigation: 
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1. Construction of bifurcation diagram and description of the topological 

structure of the system. 

2. Finding action- angle variables and explicit formulae for orbits. 

3. Quantization of the system. 

4. Painleve analysis. 

We constructed this new system using methods based on an observation 
from [9 . We will describe the methods and how we applied them elsewhere. 
It seems that it is the only new integrable system with cubic integral on 
a closed manifold containing elementary functions only that can be found 
using these methods. But still one can construct many new local examples 
of metrics with integrable geodesic flows, which we will do in a joint paper 
with Tabov and Topalov. 
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